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Recent years have seen great progress in our understanding of the electronic properties of nano¬ 
materials in which at least one dimension measures less than 100 nm. However, contacting true 
nanometer scale materials such as individual molecules or nanoparticles remains a challenge as even 
state-of-the-art nanofabrication techniques such as electron-beam lithography have a resolution of 
a few nm at best. Here we present a fabrication and measurement technique that allows high sensi¬ 
tivity and high bandwidth readout of discrete quantum states of metallic nanoparticles which does 
not require nm resolution or precision. This is achieved by coupling the nanoparticles to resonant 
electrical circuits and measurement of the phase of a reflected radio-frequency signal. This requires 
only a single tunnel contact to the nanoparticles thus simplifying device fabrication and improving 
yield and reliability. The technique is demonstrated by measurements on 2.7 nm thiol coated gold 
nanoparticles which are shown to be in excellent quantitative agreement with theory. 


Electrical characterization of metallic nanoparticles us¬ 
ing dc transport methods requires both a source and 
drain electrode to pass a current as well as a gate elec¬ 
trode to vary the chemical potential, see Fig. la. The 
most commonly used fabrication techniques are electro¬ 
migration, where a metallic wire is controllably loaded 
to failure such that a nanosize gap opens up in which 
particles might be trapped, or electrode platin g i n which 
a predefined gap is reduced electrochemically [lr[3j]. Al¬ 
though these methods are certainly feasible, the need for 
two tunnel contacts makes the fabrication complex and 
challenging and the quality and yield of devices highly 
variable. This is particularly problematic for the small¬ 
est nanoparticles - or molecules - less than around 3 nm 
in size. 

Here we present a different technique that allows the 
study of arbitrarily small nanoparticles using only a sin¬ 
gle tunnel contact. Importantly, the technique does not 
require nm positioning accuracy which facilitates the use 
of standard lithography methods. As illustrated in Fig. 
lb, we achieve this by placing the nanoparticles on a 
metal electrode and and by capacitively coupling the 
nanoparticle to a second electrode, which functions both 
as an ac reference ground and as a dc gate electrode. The 
devices are embedded in a resonant LC circuit and radio¬ 
frequency reflectometry is used as the measurement tool 
[J-Q. This approach, which can be considered as capac¬ 
itance spectroscopy |§, [g] taken into the radio-frequency 
domain, not only simplifies device fabrication but also 
offers high sensitivity and high bandwidth readout. 

The resonant circuits used in our experiments typically 
have a resonant frequency /o in the 300-500 MHz range 


which is given by the total device capacitance Cs - in¬ 
cluding parasitics - and a chip inductance L placed on the 
sample holder such that f 0 = 1/(27 t^/LGs), see Fig. Id. 
The key idea is that the capacitance depends in part on 
the ability of electrons to move on and off the nanopar¬ 
ticle at the driving frequency which in turn depends on 
the chemical potential of the nanoparticle which can be 
tuned by the gate electrode. The gate-dependent capaci¬ 
tance is probed by measuring the phase shift of a reflected 
radio-frequency signal at the resonant frequency /o- 

Before discussing experimental data, using Au 
nanoparticles, we will first examine, quantitatively, the 
expected phase response of the reflectometry technique 
for a generic ’electron-in-a-box’ system or quantum dot 
with discrete energy spectrum, coupled to a metallic lead 
at some finite temperature T, as illustrated in Fig. 2. 
In the presence of a radio-frequency signal, as in Fig. 
lb, the energy levels on the quantum dot will oscillate 
with respect to the lead. However, it is only possible 
for electrons to move on and off the quantum dot when 
an energy level is aligned with the thermally broadened 
electrochemical potential of the lead. In this case, the 
occupation probability of the level varies with the drive 
frequency and the movement of charge on and off the elec¬ 
trode in response to the rf drive amplitude can be param¬ 
eterized by an effective admittance Y = 1/ AR + joj AC 
where w = 2irf is the radial frequency and with capac¬ 
itance AC and resistance A R as shown in Fig. 2b. If 
the thermal energy is larger than the lifetime broaden¬ 
ing of the resonances, h'y, where 7 is the tunnel 

rate, it is relatively straightforward, see supplementary 
information, to obtain the following relations using rate 
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Figure 1: (a) Schematic of a conventional direct-current measurement. Two transport electrodes that are precisely aligned 
and separated by a few nanometers are required to enable electron transport through a nanoparticle, (b) Schematic of the 
radio-frequency reflectometry technique. The nanoparticle is only tunnel coupled to a single electrode. A second electrode, 
which does not need to be precisely aligned, functions both as an ac ground and dc gate electrode. The device is embedded in a 
LC resonant circuit which both simplifies device fabrication and allows for high sensitivity and high bandwidth measurements, 
(c) Representation of a gold nanoparticle coated in thiol chains used in our experiments, (d) Measured amplitude and phase 
response of a resonant circuit with resonant frequency ~ 305 MHz. 
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where the term —eaAV g takes into account the position 
of the discrete energy level with respect to the Fermi level 
of the lead using a lever arm a for conversion from gate 
voltage to energy SE3]. If the tunnel rate is compa¬ 
rable to the drive frequency, 7 ~ w, dissipation can be 
significant as previously observed in single-electron tun¬ 
neling devices EEl.' However, if tunnel rates are con¬ 
siderably faster than the drive frequency, that is 7 oj, 
the effective resistance A R diverges and the junction will 
be capacitive only. It is this limit that we will consider 
here. Using standard network analysis [HI, it can be 
shown, see supplementary information, that in this case 
the measured phase shift, for an under-coupled resonator, 
is related to the effective capacitance as: 


A-f> = —2QAC /Cy,- (3) 

where Q is the quality factor of the resonator given by the 
ratio of the resonant frequency and bandwidth A/ of the 
resonance. Using the relations 1-3 above we now have the 
tools to quantitatively investigate single-electron tunnel¬ 
ing of any device of nm dimensions in which the energy 
spectrum is discrete. 


Results 

Nanoparticle measurements. To demonstrate the 
technique experimentally, we have deposited thiol-coated 


Au nanoparticles of 2.7 nm diameter on a gold substrate. 
These particles were then covered by approximately 5.5 
nm of AI 2 O 3 deposited using atomic layer deposition, see 
Methods. On top of the oxide layer we defined a set of 
gate electrodes of varying sizes as illustrated by the scan¬ 
ning electron microscope (SEM) images in Fig. 3a. The 
results of measurements on two representative structures 
of different size in a He-3 cryostat with a base tempera¬ 
ture of ~ 400 mK are shown in Fig. 3b. The number of 
particles under the large 5x5 fim gate electrodes is not 
exactly known but, given an approximate density of Au 
nanoparticles of ~ 10-50 fi m~ 2 , see supplementary infor¬ 
mation, is likely to be large. Therefore, even though for 
individual particles the energy levels are expected to be 
well separated, the measurements on this gate electrode 
show an interleaving of the signal of many particles, re¬ 
sulting in the oscillatory pattern seen in the right panel 
image of Fig. 3b. For the measurements shown in the left 
panel of Fig. 3b, the gate overlap is significantly reduced - 
by a factor of 25 - and we observe more pronounced dips 
against a relatively flat background. For these devices 
the number of nanoparticles under the gate electrodes is 
sufficiently low for the observation of tunneling through 
individual electron states as we will examine in more de¬ 
tail. We also performed a series of control experiments, 
see supplementary information, in which no nanoparticles 
were deposited and in which we did not see any features 
in the measurements. 

Using the relations 2 and 3 obtained above for A<f> ver¬ 
sus AV^ we are able to compare the experimental data 
with theory. We obtain an excellent fit for individual 
resonances, see Fig. 3c, where the only free parameter is 
the lever arm a which sets both the width and depth of 
the dips. The width of the dips AU, oc T/a , or, more 
precisely, the full-width-half-maximum is 3.53 k B T/ea 
which is identical to what is expected in dc transport 
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Figure 2: (a) Schematic energy diagrams showing the discrete energy states of the nanoparticle and the temperature broadened 
electron distribution in the lead. The occupation probabilities of the states on the nanoparticle depend on the relative position 
of the states with respect to the electrochemical potential of the lead and varies from fully occupied (filled circles) to unoccupied 
(open circles). In response to the rf excitation, the relatively position of the levels oscillate at the applied frequency. Left: 
without a state in the thermal window around the electrochemical potential of the lead, no tunneling is allowed: all states 
are either fully occupied or empty. Right: when a level is aligned with the electrochemical potential of the lead, the average 
occupancy of the state oscillates at the driving frequency and charge thus moves back and forth between the particle and the 
lead (b) The device can be parameterized by an effective RC circuit as explained in the main text. 


of individual quantum states EMI- However, unlike 
dc transport in which the magnitude of the conductance 
is set by the tunnel rate and asymmetry of the source 
and drain tunnel junctions, in the reflectometry measure¬ 
ments - in which there is only a single tunnel junction - 
the magnitude A$ oc a 2 /T. We thus have two relations 
for a and T which can be obtained independently. In¬ 
deed, while in our experiments the temperature is known, 
T ~ 400 mK, if the temperature is set as an additional 
fit parameter, the fit procedure yields the correct experi¬ 
mental result. This enables an energy calibration without 
the need of a source-drain bias voltage and allows mea¬ 
surements of, for example, the charging energy of the 
nanoparticles or, when a magnetic field is applied, mea¬ 
surements of the g-factor or energy spectrum |8j, [9J . 

Nano-thermometry. A further important consequence 
is that our devices can be used as a primary thermometer, 
that is, a thermometer that allows temperature measure¬ 
ments without calibration by another thermometer. The 
advantage of using a thermometer based on nanoparticles 
is that it works over an exceptionally large temperature 
range. The lower limit of the thermometer is set by the 
requirement that lifetime broadening due to the coupling 
to the lead is less than thermal smearing: 2 hj < ksT. 
Since at the same time, the tunnel rates should be com¬ 
parable to or larger than the rf driving frequency, a prac¬ 
tical limit would be 7 > 100 MHz which corresponds to 
about 10 mK. The higher limit of the thermometer is set 
by the spacing between adjacent dips in measurements 
as in Fig. 3b. For a single nanoparticle this separation is 
given by the addition energy which exceeds several tens of 
meV [17] for small nanoparticles and which corresponds 
to room temperature operation fl8| . The energy-level 
separation A E of Au nanoparticles of 2.7 nm diameter is 
of order 10 meV such that for temperatures above around 
100 K a transition from the quantum (A E to 


the classical (A E <C ksT) transport regime is expected. 
In this case Eq. 2 is no longer valid, although the reso¬ 
nances will still be thermally broadened and can be mod¬ 
eled jl5l . [lG]. A further practical consideration is that the 
signal strength Ad* oc a 2 /T, and thus decreases with in¬ 
creasing temperature. This can be compensated for by 
increasing the lever arm a which is set by the gate ox¬ 
ide thickness and dielectric constant. Using AI 2 O 3 , we 
found ~ 5 nm to be the optimal oxide thickness for our 
devices: thicker oxides reduced the signal-to-noise ratio 
while thinner oxides affected device yield and resulted in 
breakthrough when the gate voltage was increased be¬ 
yond ~ 0.5 V. Nevertheless, the lever arm can be in¬ 
creased significantly by using high-dielectric oxides such 
as HfCU or TiC >2 which would improve signal-to-noise ra¬ 
tios by about two orders of magnitude. By a further 
improvement of the quality factor of the resonator, for 
example by reducing parasitic capacitances, we believe 
that room temperature operation is feasible. By embed¬ 
ding the device in a resonant circuit, as in this work, 
measurements are also fast and not limited by 1 /f noise 
as in conventional dc quantum dot thermometry |19t [20]. 

Discussion 

The experiments on Au nanoparticles described here 
are readily extended to other types of nanoparticles such 
as Pt, Fe, or Co nanoparticles, which have applications 
in a range of areas varying from catalysis to biomedicine 
[2ll . [22 ], single molecule magnets such as Mni 2 , Ci^Ni 
or TbPc 2 , or endohedral fullerenes which are of interest 
in spintronics or spin-based quantum information pro¬ 
cessing [ 23 M 25 II . As identified above, there is furthermore 
significant scope for improvement in the measurement 
sensitivity - and thus the operating temperature - of the 
nanodevices by optimization of the oxide dielectrics and 
quality factor of the resonators used. Chemical mod- 
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Figure 3: (a) Scanning electron microscope image of a nanoparticle device showing three gate electrodes with different overlap 
with the rf electrode, (b) Gate dependent phase response for a lxl jim (left) and 5x5 jum gate area (right). For the smallest 
gate areas fewer particles are being detected and individual resonances are observed, (c) Fitted curves to three individual dips 
observed in a small device. A constant background slope seen in all data, including control samples, has been subtracted. For 
the dips we obtain a = 0.0264, 0.0241, and 0.0217 from left to right, respectively. Since the widths of the dips scale, for a given 
temperature T, as AV g fwhm oc 1/a and their depths as A$ max oc a 2 the dips with lower a are both wider and shallower. 


ifLcation of the length of linker molecules, such as de- 
canedithiol used in this work, provides control over the 
tunnel rates while varying the nanoparticle solution pH 
and adsorption time |26l |, and surface modification of the 
electrodes (27i - l29| - or a combination of these - allows 
nanoparticle density control. Taken together, this opens 
up novel routes for the characterization and application 
of individual molecules and metallic nanoparticles. 

Methods 

Devices were fabricated on nominally undoped Si wafer 
substrates of resistivity > 100 Ohm cm patterned with 
Ti/Au electrodes using electron-beam lithography. Be¬ 
fore deposition of the nanoparticles, the electrodes were 
cleaned using acetone, ethanol and an oxygen plasma, 
followed by a further ethanol immersion to remove any 
surface oxides. The devices were then immersed in a 


toluene-based gold nanoparticle solution for one hour, 
followed by a rinse with toluene and ethanol. The 
gold nanoparticles have a 2.7 nm core diameter and 
were coated with octanethiol (length: 1.44 nm) and 
decanedithiol (1.82 nm) mixed monolayer. The devices 
were subsequently covered by approximately 5.5 nm of 
AI 2 O 3 using atomic layer deposition. In the final stage, 
a further set of Ti/Au gate electrodes were defined 
on the devices using electron-beam lithography. The 
devices typically showed no leakage between electrodes 
separated by the AI 2 O 3 layer up to around 0.5 V, at 
which point the oxide coating on the devices began to 
break down. 

The devices were mounted on a microstrip line coated 
printed circuit board sample holder with radio-frequency 
connections. Measurements were carried out at T ~ 400 
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mK in a Helium-3 cryostat, custom modified to allow 
reflectometry radio-frequency measurements, with a low 
noise cryogenic amplifier to improve the signal-to-noise 
ratio of the reflected signal. Demodulation is achieved 
by mixing the reflected rf signal with the reference 
signal. Both quadratures of the signal are detected, 
allowing measurements of both the amplitude and phase 
response. For the measurements shown in Fig. 3 a chip 


inductor L = 500 nH was used. This yielded a resonance 
frequency f 0 = 345 MHz and capacitance Cy, = 0.43 pF 
with quality factor Q ~ 59 which were the parameters 
used in the fit procedure for Fig. 3c. The dc gate voltage 
was applied to the samples via a bias-tee. 
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The following sections comprise the supplementary material to the main text. We provide the following: (i) An 
analysis of the RLC network used in the experiments; (ii) A derivation of the effective admittance of the nanoparticle 
devices, starting with the rate equations and (iii) The procedure used for the data fits and control experiments data 
in which no Au nanoparticles were present. 


I. RLC NETWORK 

A diagram of the circuit used in this work is shown in Fig. SI. The circuit consists of a capacitance C p which takes 
into account the parasitics of the components and wiring, an inductor L which consists of the chip inductor placed 
on the sample holder, and a resistor R m which takes into account the resistances in the wiring and the inductor - but 
could also double to allow for a matching resistor. The nanoparticle device, when modeled as a lumped component 
resistor A R and capacitor AC in parallel, as described in Section II below, allows the treatment of the entire system 
as a simple tank circuit. The impedance of the resonant circuit can then be described as: 


Ztot — Rm + icoL + A R || (AC + C p ) 


(1) 


Expanding out the terms and arranging into real and imaginary components, we arrive at the impedance of the circuit 
as a function of the ‘components’ of the device, A R and AC. 


Ztot — 


R„ 


A R 


1 + (wC s AR) 2 _ 
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1 + (ujCy,AR) 2 _ 


( 2 ) 


where C E = C p + AC. From here, it is possible to understand how a shift in the resistance or capacitance of the 
device will change the response seen on the circuit. We can simplify this equation by looking at how it behaves near 



Figure S 1: Schematic model of the RLC circuit used in the radio-frequency reflectometry measurements. The nanoparticle 
device is modeled as a lumped component resistor A R and capacitor AC in parallel which allows the treatment of the entire 
system as a simple tank circuit. 
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resonance, since this is the nominal operation range - as carried out in the resonant circuit model of Ref. [lj. Making 
the (generally good) approximation that ujARCy, » 1 and substituting uJq = 1 /(LG's): 


Ztot,res — Rm + T „„-1" lOjL - —— (3) 

ArtOs UOs 

where the imaginary component vanishes to zero exactly on resonance. Note that in our work A R is typically very 
large, which differs from the work by, e.g., Schoelkopf et al [2] where the resonant circuit is used to down-convert the 
resistance of a single-electron transistor (typically of order 10 kf2) to 50 f2. From here on, unless explicitly stated, 
the frequency is assumed to be on resonance. If we define Z r and Zi as the real and imaginary components of the 
impedance of the circuit, and Zq as the 5012 impedance of the wave guide connected to the tank circuit, then the 
reflection coefficient is given by: 


_ Z - Z 0 _ Z 2 + Z 2 - Z£ + 2 iZ 0 Zj 

_ Z + Z 0 (Zr + Z 0 ) 2 + Zf 1 ’ 

This, along with the calculated impedance from above, can be put together to show how the measured response of 
the tank circuit will vary as a function of effective dissipative A R and reactive AC components of the device, which 
can be in turn matched to the physics of the device itself using appropriate models. 

The amplitude and phase of the reflected signal are given by: 


V(z 2 + z 2 ~ z 2 y + (2 ZqZ& 

11 (Zr + Z 0 ) 2 + Zf 


(5) 


r 0 


= at an 


2 ZnZi 


( 6 ) 


The phase response is a clear indicator of under [Z(ui res ) > 5012] or over [Z(uj res ) < 5012] coupling. In the latter 
case, the phase change through resonance will be 0 < </> < it and the former will result in a phase shift of <f> = 2n. It 
is now possible to see how the magnitude and phase of the reflection coefficient vary with relevant factors. For the 
capacitance shift in phase of the reflection coefficient: 




dAC 


2Zn 


2QZ 0 


uC 2 (R 2 eff -Z 2 ) C^Reff-ZoY 


(7) 


where R e ff = R m + L/(ARCs) and Q = y Reff + Zq) is the loaded quality factor of the resonator. In the limit 
where R m —> 0 and A R —> oo this reduces to Eq. 3 in the main text. For the change in reflection magnitude due to 
a device resistance shift (still assuming ujARCy, » 1) : 
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2LZ n 
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„= Uo CvAR 2 (R eff + Z 0 ) 2 AR 

And finally, the change in magnitude due to a capacitance shift: 


( 8 ) 
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dAC 


-2Z 0 Q 2 
CsAR 


(9) 


From this, it can be seen that for a fixed device capacitance, the magnitude shift becomes far more sensitive to resistive 
changes for a device with smaller resistance. It should also be said that the change of reflection coefficient angle is 
negligible for changes of device resistance in the case where AR is large. 
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II. EFFECTIVE ADMITTANCE: RATE EQUATIONS 


To calculate the effective admittance of the nanoparticle device we start by assuming that we are in the incoherent 
regime such that the tunnel coupling 7 is sufficiently small: hj <C kt>T ■ The tunnel rates for a single electron to tunnel 
off and on a level at energy e relative to the Fermi energy of the tunnel electrode are given by |3i, ;jj: 


7 off = 7(1 - /(e)) ( 10 ) 

7 on = mfie) (U) 


where 77 = 2 indicates spin degeneracy and /(e) is a Fermi function. To simplify the analysis we will use 77 = 1 here 
which does not affect the results - apart from a temperature-dependent shift of the resonance center Q which is not 
important for the analysis below. In the presence of an rf drive, the position of the level relative to the Fermi energy 
of the tunnel electrode is given by: 


e(t) = —ea(V g + 


( 12 ) 


where we make use of the lever arm a = Cr/(Cl + Cr) where Cr and Cl are the geometric capacitances between the 
nanoparticle and the gate electrode and the nanoparticle and the rf electrode, respectively. Typically for our devices 
Cr -C Cl and as a result a< 1. The probabilities for the level to be either empty (Pp) or occupied (Pi) by a single 
charge are: 


dPp 
d t 
d^ 

d t 


loffPl - JonPo 

lonPo - loffPl 


Since Pp + Pl = 1 we arrive at the master rate equation: 


77 - + 7 P 1 = 7 /(e) 
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eaVRF iuit\ 

k B T ) 


(13) 

(14) 


(15) 


where we use a small excitation approximation, i.e. eaVRF <C k B T. Notice that the gate voltage term cannot 
be assumed small in comparison with thermal energies. Using G = e~CaV s )/( k BT) £ or } )rev i^y ; solution to this 
differential equation is a hypergeometric function Q: 


Pi 


( 1 \ v (, *7 , n. G eaVR F e iut \ 

Vl + G/ 21 V’ w’ w’l + G k B T ) 


This can be expressed as a sum: 


(16) 


Pi 


1 

1 + G 


i + E 


—27 

ujn\(n — 77 ) 

V CO / 


f G eaV RF e iut 

v 1 + G k^T 


n 


(17) 


This is the full response of a single electron revoir with a single lead and single gate to a small oscillation applied to 
the source. The terms in this sum die off rapidly with increasing n, so by taking the first term (n = 1), we can extract 
much of the physics of the system without too much complexity. 


Pi 


1 


\ % G eaVR F e^ 

1 _ *2 1 + G k B T 

UJ 


1 + G 


(18) 
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The expected charge drawn from a voltage source connected to the RF electrode at time t is given by Q R f( t) = eaP\ (t ). 
The expected current flowing from the voltage source is then: 


iRF{t ) 


dQRFjt) 

dt 


G u'y e 2 a 2 

(1 + G ) 2 cj — *7 k R T 


V RF e iut 


The effective impedance is thus: 


Z e ff = V RF e lut /i RF = 

Equating this with Z~^ = R~ff + jwC e // we find that: 


(1 + G) 2 k B T lo 2 + 7 2 
G e 2 a 2 loj( ui + ij) 


Reff ~ 

Ceff = 


ML (iL, A ( 1 + G ) 2 

e 2 a 2 7 Vw 2 ) G 
e 2 a 2 / uj 2 \ 1 G 

kLrW + ) (i + G ) 2 


Identifying (1 + G) 2 /G as 4 cosh 2 ^ 


we arrive at Eq. 1 and 2 in the main text. 


(19) 


( 20 ) 


( 21 ) 

( 22 ) 


III. DATA ANALYSIS AND CONTROL DEVICES 

Using Eq. 7 and Eq. 22 (or Eq. 2 and Eq. 3 in the main text) we are able to fit the measured data to theory, provided 
Cy. and Q are known. We obtain both from the measured amplitude and phase response over a large frequency range 
which shows a resonance at / = 345 MHz. Knowing the value of the chip inductor L = 500 nH this yields Cy = 0.43 
pF. The quality factor is measured by fitting the phase response to: = (f >o + 2 atan[2Q(l — f/fo)] where <j >o 

is the angle at the resonance frequency /q. This yielded a best fit to the data of Q = 59 ± 4. This is somewhat 
larger than expected for a Zq = 50 Q loaded circuit which we tentatively attribute to an impedance mismatch at the 
cryogenic amplifier connected to the circuit. For the data fits in the main text we used 7 u> and the experimentally 
determined Gy = 0.43 pF and Q = 59. A constant background slope, seen in all data including control samples, see 
e.g. Fig. S2b below, has been subtracted. 

The full-width-half-maximum (FWHM) of the resonances as a function of gate voltage (V g ) described by 
Eq. 22 provides a relation between temperature and lever arm: 


knT knT 

A V g ^ FWH M = -?—(ln( 3 + v^8) - ln( 3 - v^)) « 3.53-^- 


(23) 


The gate dependency is the differential of the fermi function, as expected. Using Eq. 7 and 22, the maximum phase 
signal is given by: 


|A$ 


max 


2Q e 2 a 2 
Cy 4:k R T 


From these relations we can determine the temperature and lever arm as: 


2C'e| A<I> 

max FWHM 

(3.53 ) 2 k B Q 


2U e |A$ 

max |A Vg^FWHM 

3.53eQ 


(24) 


(25) 


We are thus able to determine the temperature from measurements of the peak width and height directly which 
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Figure S 2: (a) Scanning electron micrograph of a device with Au nanoparticles as indicated by the red circles, (b) Phase 
data from a control device without nanoparticles for large 5x5 (jm gate area. The vertical scale bar (phase) is set to allow for 
a direct comparison with Fig. 3b in the main text, (c) Schematic of the nanoparticle (top) and control (bottom) devices. 


does not require a separate measurement of the lever arm. The signal strength is maximized by improving the 
lever arm and quality factor. As described in the main text, the devices can be used as a primary thermometer, 
the limits of which are set by the tunnel coupling 7 and the 0D level spacing and charging energies of the nanoparticles. 

To ensure that we measured the nanoparticles and not some other effect, e.g. charge traps in the oxides, control 
samples were prepared. When nanoparticles are deposited, a rough indication of the average density is obtained 
from scanning electron micrographs, see Fig. S2a. For the control devices, exactly the same fabrication procedure 
as outlined in the Methods section in the main text was carried out, with the exception that no nanoparticles were 
deposited onto the substrates. All other variables were kept the same, and the various control devices measured. An 
example is shown in Fig. S2b, for which indeed no signal is observed apart from a slowly decreasing background 
seen in all devices - leading to the conclusion that the nanoparticles are indeed responsible for resonances in the 
data. Samples with Au nanoparticles but with considerable thicker gate oxides (> 10 nm AI 2 O 3 ) also did not show a 
response at 400 mK (not shown) as expected since for these devices the lever arm will be much smaller and thus the 
measured signal strength, as expressed by Eq. 24. 
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